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Abstract 

Understanding the structure of the Internet graph is a crucial step for building accurate network 
models and designing efficient algorithms for Internet applications. Yet, obtaining its graph structure 
is a surprisingly difficult task, as edges cannot be explicitly queried. Instead, empirical studies rely on 
traceroutes to build what are essentially single-source, all-destinations, shortest-path trees. These trees 
only sample a fraction of the network's edges, and a recent paper by Lakhina et al. found empirically 
that the resuting sample is intrinsically biased. For instance, the observed degree distribution under 
traceroute sampling exhibits a power law even when the underlying degree distribution is Poisson. 

In this paper, we study the bias of traceroute sampling systematically, and, for a very general class of 
underlying degree distributions, calculate the likely observed distributions explicitly. To do this, we use 
a continuous-time realization of the process of exposing the BFS tree of a random graph with a given 
degree distribution, calculate the expected degree distribution of the tree, and show that it is sharply 
concentrated. As example applications of our machinery, we show how traceroute sampling finds power- 
law degree distributions in both 5-regular and Poisson-distributed random graphs. Thus, our work puts 
the observations of Lakhina et al. on a rigorous footing, and extends them to nearly arbitrary degree 
distributions. 



1 Introduction 

Owing to the great importance of the Internet as a medium for communication, a large body of recent work 
has focused on its topological properties. Perhaps most famously, Faloutsos et al. [9] exhibited a power-law 
degree distribution in the Internet graph at the router level (i.e., the level at which the Internet Protocol 
(IP) operates). Similar results were obtained in [10, 2]. Based on these and other topological studies, it is 
widely believed that the Internet's degree distribution has a power-law form with exponent 2 < a < 3, i.e., 
the fraction aj, of vertices with degree k is proportional to k~ a . These results have motivated both the search 
for natural graph growth models that give similar degree distributions (see for instance [8]) and research into 
the question of how the topology might affect the performance of Internet algorithms and mechanisms (for 
instance [17]). 
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However, unlike graphs such as the World Wide Web [13] in which links from each site can be readily 
observed, the physical connections between routers on the Internet cannot be queried directly. Without 
explicitly knowing which routers are connected, how can one obtain an accurate map of the Internet? 
Internet mapping studies typically address this issue by sampling the network's topology using traceroutes: 
packets are sent across the network in such a way that their paths are annotated with the IP addresses of the 
routers that forward them. The union of many such paths then forms a partial map of the Internet. While 
actual routing decisions involve multiple protocols and network layers, it is a common assumption that the 
packets follow shortest paths between their source and destination, and recent studies show that this is not 
far from the truth [15]. 

Most studies, including the one on which [9] is based, infer the Internet's topology from the union of 
traceroutes from a single root computer to a larger number of (or all) other computers in the network. If 
each edge has unit cost plus a small random term, the union of these shortest paths is a BFS tree. 1 This 
model of the sampling process is admittedly an idealization for several reasons. First of all, most empirical 
studies only use a subset of the valid IP addresses as destinations. Secondly, for technical reasons, some 
routers may not respond to traceroute queries. Thirdly, a single router may annotate different traceroutes 
with different IP addresses, a problem known as aliasing. These issues are known to introduce noise into the 
measured topology [1, 5]. 

However, as Lakhina et al. [14] recently pointed out, traceroute sampling has a more fundamental bias, 
one which is well-captured by the BFS idealization. Specifically, in using such a sample to represent the 
network, one tacitly assumes that the sampling process is unbiased with respect to the parameters under 
consideration, such as node degrees. However, an edge is much more likely to be visible, i.e., included 
in the BFS tree, if it is close to the root. Moreover, since in a random graph, high-degree vertices are 
more likely to be encountered early on in the BFS tree, they are sampled more accurately than low-degree 
vertices. Indeed, [14] showed empirically that for Erdos-Renyi random graphs G(n,p) [7], which have a 
Poisson degree distribution, the observed degree distribution under traceroute sampling follows a power law, 
and this has been verified analytically by Clauset and Moore [6]. In other words, the bias introduced by 
traceroute sampling can make power laws appear where none existed in the underlying graph! Even when the 
underlying graph actually does have a power-law degree distribution k~ a , Pctcrmann and De Los Rios [22] 
and Clauset and Moore [6] showed numerically that traceroute sampling can significantly underestimate its 
exponent a. 

This inherent bias in traceroute sampling (along with the fact that no alternatives are technologically 
feasible at this point) raises the following interesting question: Given the true degree distribution {a^} of the 
underlying graph, can we predict the degree distribution that will be observed after traceroute sampling? 
Or, in pure graph-theoretic terms: can we characterize the degree distribution of a BFS tree for a random 
graph with a given degree distribution? 

Our answers to these questions quantify precisely the bias introduced by traceroute sampling, while 
verifying formally the empirical observations of Lakhina et al. [14]. In addition, they can be considered as a 
significant first step toward a much more ambitious and ultimately more practical goal of inferring the true 
underlying distribution of the Internet from the biased observation. 

Our Results 

Our main result in this paper is Theorem 2, which explicitly characterizes the observed degree distribution 
as a function of the true underlying distribution, to within sharp concentration. When we say that {at} is a 
degree distribution, what we mean precisely is that the graph contains au ■ n nodes of degree k. In proving 
the result, we restrict our attention to underlying distributions which are "not too heavy-tailed," and in 
which all nodes have degree at least 3: 

Definition 1. A degree distribution {cik} is reasonable if = for k < 3 ; and there exist constants a > 2 
and C > such that < C ■ k~ a for all k. 

1 Several studies, including [2, 21], have used traceroutes from multiple sources. However, the number of sources used is quite 
small (to our knowledge, at most 12). 
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The requirement that the degree distribution be bounded by a power law k~ a with a > 2 is made mostly 
for technical convenience. Among other things, it implies that the mean degree 6 = ^2 k ka^ of the graph 
is finite (although the variance is infinite for a < 3). Note that this requirement is consistent with the 
conjectured range 2 < a < 3 for the Internet [9, 10]. The requirement that the minimum degree be at least 3 
implies, through a simple counting argument, that the graph is w.h.p. connected. 2 This is convenient since it 
ensures that the breadth-first tree reaches the entire graph. However, as we discuss below, this requirement 
can be relaxed, and in the case of disconnected graphs such as G(n,p = 6/n), we can indeed analyze the 
breadth-first tree built on the giant component. 

In order to speak precisely about a random (multi)graph with a given degree sequence, we will use the 
configuration model [3]: for each vertex of degree k, we create k copies, and then define the edges of the 
graph according to a uniformly random matching on these copies. Our main result can then be stated as 
follows: 

Theorem 2. Let {aj} be a reasonable degree sequence. Let G be a random multigraph with degree distribution 
{aj}, and assume that G is connected. Let T be a breadth-first tree on G, and let Aj bs be the number of vertices 
of degree j in T. Then, there exists a constant ( > such that with high probability, \A° hs — a° bs n| < n 1_< * 
for all j, where 



J rf- 1 m ^ Pvis (ty n (i - PvisWr 1 -™ dt 
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We can use the notion of generating functions [25] to obtain a more concise expression of all 0,^+1 as follows : 
if g(z) — Y^jLo a i z ^ > then a° bs is the coefficient of z 
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The bulk of this paper, namely Sections 2-5, is devoted to the proof of Theorem 2. In Section 6, we 
apply our general result to 5-regular graphs and graphs with Poisson degree distributions. In both cases, 
we find that the observed degree distribution follows a power law k~ a with exponent a = 1. In the case of 
Poisson degree distributions, our work thus subsumes the work of Clauset and Moore [6]. 

The proof of this result is based on a process which gradually discovers the BFS tree (see Section 2). 
By mapping it to a continuous-time process analogous to Kim's Poisson cloning model [12], we can avoid 
explicitly tracking the (rather complicated) state of the FIFO queue that arises in the process, and in 
particular the complex relationship between a vertex's degree and its position in the queue. This allows us 
to calculate the expected degree distribution to within o(l) in Section 3. In Section 5, we see how these 
calculations can be rephrased in terms of generating functions, to yield the alternate formulation of Theorem 
2. The concentration part of the result, in Section 4, analyzes a different, and much more coarse-grained, 
view of the process. By carefully conditioning on the history of the process, we can apply a small number of 
Martingale-style bounds to obtain overall concentration. 



2 A Continuous-Time Process 
2.1 Breadth-First Search 

We can think of the breadth-first tree as being built one vertex at a time by an algorithm that explores 
the graph. At each step, every vertex in the graph is labeled explored, untouched, or pending. A vertex is 

2 We say that a sequence of events £ n occurs with high probability (w.h.p.) if Prob[£ n ] = 1 — o(l) as n — > oo, and with 
overwhelmingly high probability (w.o.h.p.) if Prob[£ n ] = 1 — o(n~ c ) for all c. Note that by the union bound, the conjunction 
of a polynomial number of events, each of which occurs w.o.h.p., occurs w.o.h.p. 
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explored if both it and its neighbors arc in the tree; untouched if it is still outside the tree; and pending if 
it is on the boundary of the tree, i.e., it may still have untouched neighbors. Pending vertices are kept in 
a queue Q, so that they are explored in first-in, first-out order. The process is initialized by labeling the 
root vertex pending, and all other vertices untouched. Whenever a pending vertex is popped from Q and 
explored, all of its currently untouched neighbors are appended to Q, and the connecting edges are visible. 
On the other hand, edges to neighbors that are already in the queue are not visible. 

For the analysis in this paper, it is convenient to think of the algorithm as exploring the graph one copy 
at a time, instead of one node at a time. The queue will then contain copies instead of vertices. At each 
step, the partner v of the copy u at the head of the queue is exposed, and both of them are removed from the 
matching. Also, all of v's siblings are added to the queue, unless they were in the queue already. (We refer 
to two copies of the same vertex as siblings.) We will say that an unexposed copy is enqueued if it is in Q, 
and untouched if it is not. Thus, a copy is untouched if its vertex is, and enqueued if its vertex is pending 
and the edge incident to it has yet to be explored. Formally, the breadth-first search then looks as follows: 

Algorithm 1 Breadth-First Search at the Copy Level 
1: while Q is nonempty do 
2: Pop a copy u from the head of Q 
3: Expose u's partner v 
4: if v is untouched then 
5: Add the edge (u, v) to T 
6: Append v's siblings to Q 
7: else 

8: Remove v from Q 
9: end if 
10: end while 



An edge will be visible and included in T if, at the time one of its endpoints reaches the head of the 
queue, the other endpoint is still untouched. 

2.2 Exposure on the fly 

Because G is a uniformly random multigraph conditioned on its degree sequence, the matching on the copies 
is uniformly random. By the principle of deferred decisions [20], we can define this matching "on the fly," 
choosing -u's partner v uniformly at random from among all the unexposed copies at the time. 

One way to make this random choice is as follows. At the outset, each copy is given a real-valued index 
x chosen uniformly at random from the unit interval [0, 1]. Then, at each step, w's partner v is chosen as the 
unexposed copy with the the largest index. Thus, it is convenient to think of the algorithm as taking place 
in continuous time, where t decreases from 1 to 0: at time t, the copy at the head of the queue is matched 
with the unexposed copy of index t. Since the indices of w's siblings are uniformly random, while conditioned 
on being less than t, this approach maintains the following powerful kind of uniform randomness: at time t, 
the indices of the unexposed copies, both inside and outside the queue, are uniformly random in [0, t). 

We define the maximum index of a vertex to be the maximum of all its copies' indices. At any time i, 
the untouched vertices are precisely those whose maximum index is less than t, and the explored or pending 
vertices (whose copies are explored or enqueued) are those whose maximum index is greater than t. This 
observation allows us to carry out an explicit analysis without having to track the (rather complicated) state 
of the system as a function of time. 

At a given time t, let C uncx (t) and C unto {t) denote the number of unexposed and untouched copies, 
and let V unto ,j{t) denote the number of untouched vertices of degree j; note that C unto {t) — jV unto j(t). 
We start by calculating the expectation of these quantities. The probability that a vertex of degree j has 
maximum index less than t is exactly t 3 ; therefore, E[V^i n t ,j(i)] = ajt 3 n, and 

E[C unto (t)] = J2 3 j a 3 pn = : c unto (t)-n. (2) 
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To calculate E[C une x(£)L recall that the copy at the head of the queue has a uniformly random index 
conditioned on being less than t. Therefore, the process forms a matching on the list of indices as follows: 
take the indices in decreasing order from 1 to 0, and at time t match the index t with a randomly chosen 
index less than t. This creates a uniformly random matching on the 8n indices. Now, note that a given 
index is still remaining at time t if both it and its partner are less than t, and since the indices are uniformly 
random in [0, 1] the probability of this is t 2 . Thus, the expected number of indices remaining at time t is 

E[C uncx (t)] = 5t 2 n =: c uncx (t) ■ n . (3) 

The following lemma shows that Vuntoj^), C unt o(t) and C uncx (t) are concentrated within o(n) of their 
expectations throughout the process. Note that we assume here that the graph G is connected, since 
otherwise, the process is not well-defined for all t e [0, 1]. 

Lemma 3. Let {cij} be a reasonable degree distribution, and assume that G is connected. Then, for any 
constants (3 < min(i, and e > 0, the following hold simultaneously for all t 6 [0, 1] and for all j < n, 

w.o.h.p.: 

\V unto ,j{t) - aj t j -n\ < n 1/2+£ 
|C uncx (t) - c uncx (t) • n\ < n 1/2+t 

|Cunto(i) - Cunto(i) ' «| < Tl 1 ' ~ l \ 

where c unto {t) and c UDCX (t) are given by (2), (3). 

Note that this concentration becomes weaker as a — > 2, since then [3 — > 0. 
Proof. Our proof is based on the following form of the Hoeffding Bound [11, 16]: 

Theorem 4 (Theorem 3 from [16]). If X\, . . . ,Xk are independent, non-negative random variables with 
Xi < bi for all i, and X = J2i %i> then for any A > 0: 

Prob[|X-E[X]| > A] < 2c- 2A2 /£» fc - . 

First, V un to,j{t) is a binomial random variable distributed as Bin^n, V). By applying Theorem 4 to aj-n 
variables bounded by 1, the probability that V un to,j(t) differs by A = n 1,/2+e from its expectation is at most 
2e~ 2n2e / aj < e - ™ 2 \ Thus, at each individual time t and for each j, the stated bound on V unto j(t) holds 
w.o.h.p. 

We wish to show that this bound holds w.o.h.p. for all j and all t, i.e., that the probability that it is 
violated for any t and any j is o(n~ c ) for all c. Notice that the space of all times t is infinite, so we cannot 
take a simple union bound. Instead, we divide the interval [0, 1] into sufficiently small discrete subintervals, 
and take a union bound of those. Let m = • jajn — Sn be the total number of copies, where 6 is the mean 
degree (recall that 8 is finite, because {aj} is reasonable). We divide the unit interval [0, 1] into m h intervals 
of size m~ b , where b will be set below. By a union bound over the (™) pairs of copies, with probability at 
least 1 — to 2-6 , each interval contains the index of at most one copy, and therefore at most one event of 
the queue process. Conditioning on this event, V un to,j(t) changes by at most 1 during each interval, so if 
Kintoj(i) is close to its expectation at the boundaries of each interval, it is close to its expectation for all 
t E [0, 1]. In addition, we take a union bound over all j. The probability that the stated bound is violated 
for any j in any interval is then at most 

n(m b e-" 2e +m 2 - 6 ) = 0(n 3 - b ) , 

which is o(n~ c ) if b > c + 3. 

For the concentration of C unex (i), we notice that unexposed copies come in matched pairs, both of 
which have index less than t. Therefore, C uncx (t) is twice a binomial random variable distributed as 
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Bin(J^ jajn/2, t 2 ). Applying Theorem 4 with A = n 1//2+£ gives the result for fixed t, and taking a union 
bound over t as in the previous paragraph shows the concentration of C unox (t). 

To prove concentration of C un to(i) for fixed t, we let X, be the number of copies of node i that are 
untouched at time t. Then, C unto (t) — J2i-%i, and the denominator in the exponent for the bound of 
Theorem 4 is 

E 6 ? = J2f a i n < CnJ2f' a < { O(nlogn) a = 3 



0{n) a > 3 



Hence, by Theorem 4, whenever (3 < min(|, ^y 2 -), we obtain that |C unto (t) - E[C un to(*)]| < n 1 13 w.o.h.p. 
A union bound over t as before completes the proof. ■ 



3 Expected degree distribution 

In this section, we begin the proof of Theorem 2 by analyzing the continuous-time process defined in Section 2, 
and calculating the expected degree distribution of the tree T. 

By linearity of expectation, the expected number of vertices of degree j in T is the sum, over all vertices 
v, of the probability that j of v's edges are visible. Consider a given vertex v of degree i. It is touched when 
its copy with maximum index is matched to the head of the queue, at which time its i — 1 other copies join 
the tail of the queue. If to of these give rise to visible edges, then v's degree in T will bera+1, namely these 
to outgoing edges plus the edge connecting v back toward the root of the tree. 

Let pi^m denote the probability of this event, i.e., that a vertex of degree i has to copies that give rise to 
visible edges. Then the expected degree distribution is given by 

i 

Moreover, let Pi, m (t) denote the probability of this event given that v has maximum index t. Then, since 
t is the maximum of i independent uniform variables in [0, 1], its probability distribution is di l /dt = it' -1 , 
and we have 

Pi,m = / if' 1 p it m(t) dt . (5) 
Jo 

Our goal is then to calculate /Oj lTO (t). 

Let us start by calculating the probability P v i s (t) that, if v has index t, a given copy of v other than the 
copy with index t — that is, a given copy which is added to the queue at time t — gives rise to a visible edge. 
Call this copy u, and call its partner w. According to Algorithm 1, the edge (u, w) is visible if and only if (1) 
u makes it to the head of the queue without being matched first, and (2) when it does, w is still untouched. 
But (1) is equivalent to saying that w is untouched at time t, since if to is already in the queue at time t, it 
is ahead of u, and u will be matched before it reaches the head of the queue. Similarly, (2) is equivalent to 
saying that all of w 's siblings' partners are untouched at time t, since if any of these are already in the queue 
at time t, and thus ahead of u, then to's vertex will be touched, and w enqueued, by the time u reaches the 
head of the queue. 

Given the number of untouched and unexposed copies C un to(t) and C unex (t) at the time t when u joins 
the queue, the probability that its uniformly random partner w is untouched is P un to(t) = C unto (t) / C unex (t) . 
Conditioning on this event, the probability that w belongs to a vertex with degree k is P U nto,fe(*) = 
kVunto,k(t)/C unto (t). We require that the partners of to's k — 1 siblings are also untouched. If we ignore the 
fact that we are choosing untouched copies without replacement (and that one untouched copy has already 
been taken for to), and if we assume that v, its neighbors, and its neighbors' neighbors form a tree (i.e., that 
v does not occur in a triangle or 4-cycle, and that neither it nor its neighbors have any multiple edges), then 
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the probability that these k — 1 copies are all untouched is P un to(t) k 1 - This gives 

Pvis(t) = Punto(t) Punto,k(t) P\into(t) k 1 

k 

= ^Punto,fc(t)Punto(*) fe . (6) 
k 

Since Kmto,fc(i)i C un t (t) and C uncx (t) are concentrated according to Lemma 3, substituting their expec- 
tations then gives a good approximation for P v ; s (t), namely 

p (t) = kaktk ( Cunt °ft) \ 

, c unto(^) \ c uncx(^)/ 

Then, if we neglect the possibility of self- loops and parallel edges involving u and its siblings, and again 
ignore the fact that we are choosing without replacement (i.e., that processing each sibling changes C unto , 
Cunex, and Pvis slightly) the events that each of u's siblings give rise to a visible edge are independent, and 
the number m of visible edges is approximately binomially distributed as Bm(i — l,p v i s (t)). 

We wish to confirm this analysis by showing that w.h.p. u, its neighbors, and its neighbors' neighbors form 
a tree. It is easy to show this for graphs with bounded degree; however, for power-law degree distributions 
ak ~ k~ a , it is somewhat delicate, especially for a close to 2. The following lemmas show that there are 
very few vertices of very high degree, and then show that the above is w.h.p. true of u if u has sufficiently 
low degree. We then show that we can think of all the copies involved as chosen with replacement. Recall 
that the mean degree S — J2j 3 a j i s finite, and let < min(i, ^ a ~ 2 ^ ) as in Lemma 3. 

Lemma 5. The probability that a random copy belongs to a vertex of degree greater than k is o{k~ 2 P). 

Proof. This probability is 

^> fcJQj < 9. V i 1 '" < - h- (a - 2) - o(k- 213 ) 
Zjjai < S(2-a) k ~ « k 



Lemma 6. There are constants 7 > r\ > such that if v is a vertex of degree i < n v , then the probability 
that v or its neighbors have a self-loop or multiple edge, or that v is part of a triangle or a cycle of length 4 7 
is o(n -7 ). Thus, v, its neighbors, and its neighbors' neighbors form a tree with probability 1 — o(n -7 ). 

Proof. First, we employ Lemma 5 to condition on the event that none of u's neighbors have degree greater 
than n x , where A (and rf) will be determined below. By a union bound over these i < n v neighbors, this 
holds with probability 1 — o(n J) ~ 2A/3 ). (Unfortunately, we cannot also condition on u's neighbors' neighbors 
having degree at most n x without breaking this union bound.) 

Now, if we choose two copies independently and uniformly at random, the probability that they are both 
copies of a given vertex of degree j < n x is j(j — l)/(Sn) 2 < n 2X ~ 2 , and the probability that they are both 
copies of any such vertex is at most n 2X ~ 1 . Moreover, the probability that two random copies are siblings, 
regardless of the degree of their vertex, is 



sib 



Taking a union bound over all pairs of copies of v, the probability that u has a multiple edge, i.e., that two 
of its copies are matched to copies of the same neighboring vertex, is at most z 2 n 2A_1 = 0(n 2v+2X ~ 1 ), and 
the probability that u contains a self-loop, i.e., that two of its copies are matched, is 0{i 2 /(Sn)) = 0(n 2l ' _1 ). 
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For each of v's neighbors, the probability of parallel edges involving it is at most n 2A P s ib = o(n 2A_2/3 ), and 
the probability of a self-loop is 0(n 2X / (5n)) = <3(n 2A_1 ). Taking a union bound over all of v's neighbors, 
the probability that any of them have a self-loop or multiple edge is o(n n+2X ~ 2/3 ). 

To determine the expected number of triangles containing v, we notice that any such triangle contains 
two copies each from v and two of its neighbors, and edges between the appropriate pairs. A given pair of 
copies is connected with probability 0(l/(Sn)), so the expected number is 

0{n 2 n^(n 2X ) 2 /(5nf)=0{n^ +AX - 1 ) . 

Similarly, each 4-cycle involves two copies each of v and two of its neighbors, such that one copy from 
each of the neighbors is matched with one copy of v, and the other two copies are matched with copies of 
the same node. Thus, the expected number of 4-cycles involving v is 

0(n 2 n 2f '(n 2A ) 2 P sib /(<5n) 2 ) =o(n 2 "+ 4A - 2 ' 3 ) . 

Collecting all these events, the probability that the statement of the lemma is violated is 

o(n~ 7 ) where 7 = — max(?7 — 2X/3, 2i] + 4A — 2(3) . 

If we set rj = (3 2 /6 and A = (3/4, then 7 = (3 2 /3. ■ 

The next lemma shows that, conditioning on the event of Lemma 6, the copies discussed in our analysis 
above can be thought of as chosen with replacement, as long as we are not too close to the end of the process 
where untouched copies become rare. Therefore, the number of visible edges is binomially distributed. 

Lemma 7. Let ij,-f be defined as in Lemma 6. There exists a constant 9 > such that for t E [n~ e , 1] and 
i < n-i, 

\pi, m {t) - Prob[Bin(i - l,P vis (i)) = m]| < n" 7 . 

Proof. Let d m ; n be the minimum degree of the graph, i.e., the smallest j such that aj > 0. Note that 
dmin > 3, and set = (3/(2d min ) < 1/12. For t > , we have that E[C uncx (t)] = 5t 2 n = fXn 1 -^"""), and 
this bound holds w.o.h.p. by Lemma 3. 

Conditioning on v's neighbors having degree at most n A as in Lemma 6, the number of visible edges of v 
is determined by a total of at most n v+x copies. These are chosen without replacement from the unexposed 
copies. If we instead choose them with replacement, the probability of a collision in which some copy is 
chosen twice is at most (n^ +x ) 2 /C uncx (i) = 0( n 2^+2A+ / 3/d ml „-i- ) = ( n -i/2) xhis can be absorbed into the 
probability o(n~ 7 ) that the statement of Lemma 6 docs not hold. If there arc no collisions, then we can 
assume the copies are chosen with replacement, and each of v's i — 1 outgoing edges is independently visible 
with probability P v i s (t), as defined in (6). ■ 

The next three lemmas then show that p v i s (t) is a very good approximation for P v i s (t) for most t, and 
that therefore the distribution of Pi, m (t) is very close to Bin(i — l,p v ; s (t)). 

Lemma 8. Let 7 and 6 be defined as in Lemma 6 and Lemma 7. There exists a constant k > such that 
for allte [n~ 9 , 1 - n~ K ], w.o.h.p. \P v i s (t) - p v is{t)\ < n~~< . 

Proof. Recall that 6 = (3/(2d m [ n ). From Lemma 3, since t > n~ e we have W.O.h.p. Cunex(^) — Q(t n) — 
f2(n 1- ' 3 / dmin ) as in the previous lemma, and C unto (i) = £l(i dmin n) = f^n 1- ' 3 / 2 ). For definiteness, take 
e = 1/12 in Lemma 3; then w.o.h.p. 

C UDCX (t) = c uncx (t)n + o{n 7 ' 12 ) = c uncx (t)n ■ (1 + o(n^ d ^- 5 / 12 )) 
C un to(0 = c unto (t)n + oin 1 -?) = c unto (t)n ■ (1 + o(n" ,3/2 )) 
Kmto,fc(*) = a k t k n + o{n 7 / 12 ) . 
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Recall that /3 < 1/2 and d min > 3. Since 0/d min - 5/12 < -1/4 < -(3/2, 

r> l.\ C un to(t) C un to(i) / , -/3/2\\ 

Now, we compare P v i s (t) with p v is(0 term by term, and separate their respective sums into the terms with 
3 < k < to 9 / 12 and those with k > n ' 12 . For all k < to 3 / 12 , since 0/12 + 0/2 - 5/12 < -1/8 < -0/4, 



kV m to,k{t) _ ka k t k | ( n -<3/4) 

Cunto(t) C unto (t) 



-^unto,A:(^) 

and since (1 + x) v = 1 + O(xy) if xy < 1, 



Cunto (^) 
Cunex(^) 



(l + o{n-^' 12 )) . 



Thus, each term obeys 



unto.fc 



(t)P mto (t) k = (^HttT + o( „-^) 



and the total error from the first n' 3 / 12 terms is at most n* 3 / 12 • o(n ^/ 4 ) = o(n ^Z 6 ). 
On the other hand, if k > n 13 ^ 12 , then for any t < 1 — n~ K we have 

i* < e"*""" < e-"" 712 "" , 

Setting k < 0/12 makes this exponentially small. In that case, taking a union bound over all rfil 12 < k < n, 
w.o.h.p. there are no unexposed vertices of degree greater than rfil 12 ; thus P un to.k(t) = and these terms of 
Pvis(i) are zero. The corresponding terms of p v is(i) are exponentially small as well, so the total error from 
these terms is exponentially small. Thus, the total error is o(n~ /3//6 ), and since 7 = 2 /3 < 0/6, this can be 
absorbed into the probability o(n~ 7 ) that the conditioning of Lemma 6 is violated. ■ 

Lemma 9. For any s,m,p and A, 

|Prob[Bin(s,p) = m] — Prob[Bin(s,p + A) = m] \ < sA . 

Proof. It is sufficient to bound the derivative of these probabilities with respect to p as follows. 

dp (to) ^ 

)p m (i- P y 



d 

— Prob[Bin(s,p) = to] 



\i- P y 

TO S — TO 



P 1-p 



< ' W -*)'-•» max ™ 

\m ) \ p 1 — p 

<™Jj2( s )-p m {i-py- m , i2( s Y—^P m (i- P y 

= max(s, s) = s . 
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Lemma 10. Let {a,j} be a reasonable degree distribution and assume that G is connected. There are constants 
0, k, rj,^i> 0, such that for all t e [n~ e , 1 — n~ K ] and all i < n v , for sufficiently large n, 

\pi, m {t) - Prob[Bin(i - l,p vis (t)) = m]\ < n~» , 

where p v i s (t) is defined in (7). 

Proof. Given Lemma 8, we apply Lemma 9 and the triangle inequality. In this case, we have s < n v and 
A < n~ 7 , so the error in pj jTO is at most rf 1 ^ 1 . Recalling from the proof of Lemma 6 that r\ = /3 2 /6 and 
7 = /3 2 /3, for sufficiently large n this is less than for any \i < (3 2 /Q. ■ 

Finally, combining Lemma 10 with (4), (5), and (7), if 



obs 
l m+l 



J2 a 

i 

where, combining (7) with (2) and (3) 



/W _1 

Jo V m 



PviB(*r(i-Pvi B (t)) < - i - m dt 



y.,.i",. ij \ 



(8) 



then wc have the following lemma. 

Lemma 11. Let {ai} be a reasonable degree sequence and assume that G is connected. There is a constant 
C > such that for sufficiently large n, for all j < n 

|E[A° bs ] -a° bs n| <n x - c . 

Proof. There are three sources of error in our estimate of E[A° bs ] for each j. These are the error n _A1 in 
Pi.m(t) given by Lemma 10, and the fact that two types of vertices are not covered by that lemma: those 
with degree greater than n v , and those which join the queue at some time t £ {n~ e , 1 — n~ K ]. The total error 
is then at most n 1 ^ plus the number of vertices of either of these types. The number of vertices of degree 
greater than n v is at most 



n 

j>n^ j>n r i 

The number of vertices that join the queue at a time t [n -8 , 1 — n~ K ] is at most the number of copies whose 
index is outside this interval. This is binomially distributed with mean n 1 ^ 6 + n 1_K , and by the Chernoff 
bound, this is w.o.h.p. less than n 1 ^^ for sufficiently large n for any £ < min(#, k). The (exponentially small) 
probability that this bound is violated can be absorbed into n 1_< » as well. Setting ( < min(/x, (a — l)r], 9, n) 
completes the proof. ■ 

4 Concentration 

In this section, we prove that the number A^ bs of nodes of observed degree j is tightly concentrated around 
its expectation E [A° bs ] . Specifically, we prove 

Theorem 12. There is a constant p > such that, with overwhelmingly high probability, the following holds 
simultaneously for all j: 

\A° hs -E[A° hs ] \ ^Oin 1 -") . 
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Proof. In order to prove concentration, the style of analysis in the previous section will not be sufficient. 
Intuitively, the reason is that changing a single edge in the graph can have a dramatic impact on the resulting 
BFS tree, and thus on the observed degree of a large number of vertices. As a result, it seems unlikely that 
A° hs can be decomposed into a large number of small contributions such that their sum can easily be shown 
to be concentrated. In particular, this rules out the direct application both of Chernoff-style bounds and of 
martingale-based inequalities. 

There is, however, a sense in which martingale bounds will prove helpful. The key is to decompose 
the evolution of A° hs into a small number of "bulk moves," and prove concentration for each one of them. 
Concretely, assume that the BFS tree has already been exposed up to a certain distance r from the root, and 
that we know the number of copies in the queue, as well as the number of untouched copies at that point. 
Since all these copies will be matched uniformly at random, one can use an edge-switching martingale bound 
to prove that the degree distributions of nodes at distance r + 1 from the root will be sharply concentrated. 
In fact, this concentration argument applies to the observed degrees of the neighbors of any "batch" of copies 
that comprise the queue Q(t) at some time t. 

We will implicitly divide the copies in the graph into such "batches" by specifying a set of a priori fixed 
points in time at which we examine the system. That is, we will approximate A° hs by the sum of the observed 
degrees of the neighbors of Q(t) over these time steps. We will show that each of the terms in the sum is 
sharply concentrated around its expectation, and then prove that the true expectation of A° hs is not very 
far from the expectation of the sum that we consider. For the latter part, it is crucial that most vertices be 
counted exactly once in the sum; this will follow readily from the concentration given by Lemma 3. 

To make the above outline precise, we let Q(t) := \Q(t)\ be the number of copies in the queue at time t, 
and let 

q(t) := E[Q(t)} = (c uncx (t) - c unto {t)) ■ n 

be its expected size. We define a sequence of r < log 2 n times at which we observe the queue and its 
neighbors. We start with t\ — 1. For each i, we let t i+1 > be maximal such that 

(ti+i) > q(ti)/n + 2n , 

where /3 is defined as in Lemma 3. Depending (deterministically) on the properties of the real- valued functions 
c uncx and c un to, there may be an i < log n such that ti+i does not exist, namely when c uncx (ti) < 2n-". If 
so, we let r be that i; otherwise, we let r = log 2 n. 

For each degree j, let Bj(i) denote the number of vertices adjacent to Q(U) whose observed degree is j. 
Lemma 13 below shows that each Bj(i) is sharply concentrated. However, we want to prove concentration 
for the overall quantity A° hs . Using a union bound over alH = 1, . . . , r and summing up the corresponding 
Bj(i) will give us concentration for A° hs , assuming that (1) not too many times U are considered, (2) nodes 
are not double-counted for multiple i, and (3) almost all nodes are considered in some batch i. 

For the first point, recall that we explicitly chose r — 0(log 2 n). For the second, observe that whenever 

C uncK (U) - C uacK (t l+1 ) > q{U) + 2n 1 -P > Q{U) 

for all times i, then all of the Q(U) are disjoint. Each of these bounds holds w.o.h.p. by Lemma 3, and by 
the union bound, w.o.h.p. they hold simultaneously. 

This leaves the third point. Here, we first bound the number of nodes that remain unexposed after time 
t r . If the construction terminated prematurely (i.e., r < log 2 n), then the fact that c uncx (0) = implies 
that c U nto(tr) < 2n -/3 , so by Lemma 3, at most 0(n 1_/3 ) copies remain unexposed w.o.h.p. On the other 
hand, when r = log 2 n, we can use the fact that the diameter of a random graph is bounded by log 2 n with 
probability at least 1 — n~ 1//2 , which we prove in the full paper using techniques of Bollobas and Chung [4]. 
Even if C U nto(tr) were fl(n) in the remaining case, since this occurs with probability at most n~ x / 2 , we have 
Cunto(tr) • n = Oin 1 / 2 ) = Oin 1 -! 3 ) since < 1/2. 

Let £ denote the event that V UYl t .j{U) — a,jt\n < n 1 / 2+e and \Q(U) — q(U)\ < 2n 1 ~P hold simultane- 
ously for all i. By Lemma 3, £ occurs w.o.h.p. In that case, we know that (1) all of the sets Q(tj) are 



11 



disjoint, and (2) the union of all the Q{U) excludes at most 2r • n 1 @ + 0(11} ^) = (^(n 1 ^) copies total 
(where O includes polylog(n) factors). Thus, \A° hs — 5Zi=i -^i (*) | = 0(n 1_/3 ) w.o.h.p., which implies that 



\e[a 



obsl 



Y^i=i E[_Bj- (i)] J = 0(n l since this difference is deterministically bounded above by n. 



By Lemma 13 below and a union bound over all i, there is a r > such that w.o.h.p. \Bj(i) — ~E[Bj(i)]\ = 
0(n 1 ~ T ) holds simultaneously for all i = 1, . . . , r and all j. Hence, by a union bound with the event £ , and 
the triangle inequality, the following holds w.o.h.p.: 



>tobs 



-E[A 



< 



i=l 



Y,\Bj(i) -mm + 



i=l 



i=l 



= 0(n 1 - p ) . 

for any p < min(/3, r), completing the proof of Theorem 12. 

The concentration for one "batch" of nodes at time t{ is captured by the following lemma. 

Lemma 13. There is a constant r > such that, for any fixed i, w.o.h.p., \Bj(i) — ~E[Bj(i) 
holds simultaneously for all j. 



= 0(n 1 - T ) 



Proof. As explained above, the idea for the proof is to apply an edge-exposure Martingale-style argument 
to the nodes that are adjacent to Q(ti). We use the following concentration inequality for random variables 
on matchings due to Wormald [26, Theorem 2.19]. A switching consists of replacing two edges {pi,P2}, 
{P3,P4> by {pi,p 3 }, {p2,Pi\- 

Theorem 14. [26] Let Xk be a random variable defined on uniformly random configurations M, M' of k 
copies, such that, whenever M and M' differ by only one switching, 

\X k (M) - X k (M')\ < c 

for some constant c. Then, for any r > 0, 

Prob[|A fe - E[X k }\ > A] < 2e- A2 /( fec2 > . 

For fixed values q and b = 61, . . . , b n , let £ q ^ denote the event that Q(U) = q and V un t j(ti) = bj for all j. 
Conditioned on £ q ,b, the matching on the q+^2j jbj copies is uniformly random. Since any switching changes 
the value of Bj(i) by at most 2, Theorem 14 implies that 



\B ] (i)-V[B 3 {i)\£ q , b }\<n 



1/2+e 



(9) 



holds w.o.h.p. for any e > 0. If we knew the queue size q and the number bj of untouched nodes of degree j 
exactly, then we could apply Theorem 14 directly. 

In reality, we will certainly not know the precise values of q and b. Therefore, we need to analyze the 
effect that deviations of these quantities will have on our tail bounds. We do this by showing in Lemma 15 
below that the conditional expectations ~E\Bj(i) \ £q t b] are close to the actual expectations F>[Bj(i)]. It follows 
that concentration around the conditional expectation implies concentration around the actual expectation. 
Specifically, write 

for the interval of possible queue lengths under consideration, and, for some < e < 1/2, write 

a^n-n^+^ajtln + n 1 ^' 

for the interval of possible numbers of untouched vertices of degree j, as well as I b := l\ x • • • x for the 
range of all possible combinations of numbers of untouched vertices. Now, let £< be the event that Q(U) € I q 
and V un to,j{ti) € Ij for all j. Notice that £< occurs w.o.h.p. by Lemma 3. 
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Lemma 15 ensures that whenever q G I q and b G I b , then the conditional expectation is close to the true 
expectation, i.e., for some r > 0, 

\E[B j (i)\£ q , b ]-E[B j (i)]\=0{n 1 - r ) . 

Thus, for all such q and b, combining this with (9) and the triangle inequality gives \Bj(i) — E[Bj(i)]\ = 
0(n 1 ~ T ), so the latter occurs w.o.h.p. Finally, a union bound with the event £< and over all j completes 
the proof. ■ 

The final missing step is a bound relating the conditional expectation of Bj(i) with its true expectation. 
Intuitively, since all relevant parameters are sharply concentrated, one would expect that the conditional 
expectation for any of the likely values is close to the true expectation. Making this notion precise turns out 
to be surprisingly cumbersome. 

Lemma 15. There is a constant r > such that, for any q G I q and b G I b , we have 

\E[B j (j)\E g , h ]--E[B j (j)]\=0(n 1 -T) . 

Proof. We first compare the conditional expectations for two "scenarios" of queue lengths and untouched 
vertices when the scenarios are close. We will see that the conditional expectations in those two scenarios will 
be close; from that, we can then conclude that any conditional expectation is close to the true expectation. 
Given q, q' and b, b', such that \q — q'\ < 4n 1 ~ /3 , and \bj — bj | < 2n 1 / 2+e for each j, we let q = min(g, q') 

and bj = mm(bj,b'j), and define the events £ := £ q .h,£' '■= £ q ' ,w , and £ := £.^. Now, we claim that, for 
some t > 0, 

E[Knto,^) I £} - v[VuntoAu) \e]\ = oin 1 --) 

for all j, and similarly for £' . By the triangle inequality, this immediately implies that 

|E[Knto,ife) I £} - E[Knto,,(^) | £']\ = Oin 1 --) . 

To prove the claim, imagine that in the (q, b) instance, we color an arbitrary, but fixed, set of q — q 
of copies in the queue black, as well as the copies of an arbitrary set of bj ~ bj vertices for each degree j. 
To expose the matching, we first expose all the neighbors of black copies, and color them blue, and then 
choose a uniform matching among the remaining (white, say) copies. The number of blue copies obeys some 
distribution -D g ,b, but in any case, it never exceeds the total number of black copies. Since q, q' G I q and 
b, b' G I b , for any v > 0, this total number is at most 

(q -?) + £ J ■ (bj -bj)< ^- f3 + 2nV^ £ j + 2 £ j . ttj n 

= An 1 -! 3 + 0(n^ 2+e + 2 ") + Oin 1 -^- 2 ^) 
= 0(n 1 " T ) , 

for any r < min(/3, 1/2 — e — 2v, (a — 2)v). Note that r > as long as 1/2 — e — 2v > 0; recall that we took 
e < 1/2 in the previous lemma, so we can choose any v < (1/2 — e)/2. 

Now, in the (q, b) instance, we can generate a uniformly random matching as follows: we choose a number 
k according to -D 9 .b, choose k copies uniformly at random and color them blue, and determine a uniformly 
random matching among the blue copies only. Then, we match up the remaining white copies uniformly at 
random. We will call a node black if at least one of its copies is black, blue if at least one of its copies is 
blue, and white otherwise. 

Since the set of nodes that are not black is deterministically the same in both instances, and the probability 
distribution of blue nodes is the same in both, the expected number of white nodes that end up with visible 
degree j is the same in both experiments. Hence, the expected total number of nodes with observed degree 
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j can only differ by the number of blue or black nodes. Even if the degrees of those nodes were chosen 
adversarially, the difference cannot be more than 0(n 1_r ), since this is a deterministic upper bound on the 
number of black or blue copies, and hence on the number of black or blue nodes. By summing up over the 
entire probability space, this now proves the claim for £ and £, and thus also for £ and £' . 

We know that if q,q' £ I q and b,b' 6 I b , then they always satisfy the necessary conditions, and hence 
the conditional expectations are within 0(n 1 ~ T ). Summing up over all q G I g and be/' therefore shows 
that 

\E[B j (i)\£ qih ]-E[B j {i)\£<]\=0(n 1 - T ) . 

Finally, because £< occurs with overwhelmingly high probability, and Bj(i) is bounded by n, we obtain that, 
for all c, 

\E[Bj(i) | £<} - E[Bj (i)] | =0(n" c ) , 
and the triangle inequality completes the proof. ■ 

5 Generating functions 

In this section, we use the formalism of generating functions [25] to express the results of Section 3 more 
succinctly, and complete the proof of Theorem 2. Given the generating function of the degree sequence of 
the underlying graph 

g{z) = ^a lZ l , 

i 

our goal is to obtain the generating function for the expected degree sequence of the breadth-first tree as 
approximated by Lemma 11, 

„ohs( „a \ A ^obs A 

9 {z)=2-^ a T- z ■ 

i 

Using the generating function formalism, we can write 

CuntoM - tg'(t), S - g'(l), c unox (t) - tV(l) , 

and from (7) we have 
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Then, combining (7) and (8), the generating function for the observed degree sequence is given by 



i-l 



i 

= z E° 



m=0 
1 



it 1 



i - 1 

m 



i - 1 
m 



m— C 

= z Vo, / it'" 1 (1 - (1 - z) Pvis (t)) i - 1 d 
= z / V oi* • [t (1 - (1 - z)Pvi S (t))] i_1 dt 

jo : 



Pvis (t) m (i- Pvis (t)) i - 1 - m di 
(^ vis (i)) ro (i-p vis (t))^ 1 - m di 



z / ff '[i(l-(l-z)p vis (i))] dt 

1 ,r i-z ,/«,'(*) 

s r?(i)n?(i) 

which completes the proof of Theorem 2. 



dt . 



Our definition of "reasonable" degree sequences implies that the graph is w.h.p. connected, so that every 
copy is eventually added to the queue. For other degree sequences, Molloy and Reed [18, 19] established 
that w.h.p. there is a unique giant component if Y]j cij(j 2 — 2j) > 0, and calculated its size within o{n). We 
omit the details, but g ohs (z) is then given by an integral from to to 1, where to is the time at which the 
giant component has w.h.p. been completely exposed; this is the time at which c unto (i) = c unox (i), namely 
the largest root less than 1 of the equation 



(11) 



6 Examples 

6.1 Regular graphs 

Random regular graphs present a particularly attractive application of the machinery developed here, as the 
generating function for a 5-regular degree sequence is simply g(z) = z s . From (1), we derive the generating 
function for the observed degree sequence: 



g obs {z) = zS ■ [ 1*- l (l-(l-z)1*V- 2 ')) s - l dt . (12) 
Jo 

This integral can be expressed in terms of the hypergeometric function 2F1 [23]. In general, for all a > — 1 
and b > 0, we have 

f\an J>^-c^ 1 T,( a + l a + b+l \ 
I t{l - Xt) At =—l* F \—' C; X ) ' 

where 

T(s + i) T(t + i) T(u) z l 



2 F 1 (s,t;u;z) = ^ 



(| T(s) T(t) T(u + i)i\ ' 
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and r(s + i)/T(s) is the rising factorial (s), = s(s + l)(s + 2) . . . (s + i — 1), also known as the Pochhammer 
symbol. In (12), a = 5 — 1, b = 5(5 — 2), and c = 1 — 5 (note a > — 1 and b > since 5 > 2) giving 

5° bs M = * • 2*1 1 - <J; 1 + 1 - . (13) 

Another useful identity is that for any negative integer q, 

2 F 1 (p, g; r; x) - ^ ^ ~ P ~ 9 j 2 F 1 (p, q; p + q + 1 - r; 1 - x) . 
Y(r -p)T(r - q) 

Here, g — c — 1 — 5, and <5 is an integer greater than 2. Thus, (13) becomes 

„ob.^_„ m + s^ms) i 



r(*) v=i i \ ^' 



r((5+ ^)( ( 5-2)^ o V -5-27 m! ' 

where the summation now ranges over only the non-zero coefficients, i.e., for all m > 5, the rising factorial 
term (1 — 5) m = and the corresponding coefficients are zero. Thus, the expected observed degree sequence 
is given by 

obs rffl r(m + ^) 

m+1 T(6+^- 2 )(5-2)m\ ■ 
To explore the asymptotic behavior of a^+i' note that 

r(m) < r(m + e) < T(m) m e 
for all m > 2 and all < e < 1. Therefore, for m > 2, we can bound o^+i as follows: 

m- 1 obs TO -i+i/(*-2) 

5V(«-2) (6 - 2) K 0m+1 < 5-2 

For any fixed 5, this gives a power-law degree sequence, and in the limit of large 5, one observes a„+i ~ m~ x ■ 
Thus, even regular graphs appear to have a power-law degree distribution (with exponent a — > 1 in the limit 
5 — > oo ) under traceroute sampling! 



6.2 Poisson degree distributions 

Clauset and Moore [6] used the method of differential equations to show that a breadth-first tree in the giant 
component of G(n,p = 5/n) has a power-law degree distribution, a m+ i ~ mT 1 for m < 5. Here, we recover 
this result as a special case of our analysis. Recall that w.h.p. the degree distribution of G(n,p — 5/n) is 
within o(l) of a Poisson degree sequence with mean 5. The generating function is then g(z) = e~ s ^~ z \ and 
the generating function for the observed degree sequence is 

g° hs (z) = z5 ■ f e-^e-^-^^^'dt 

J to 

= z I''** e-^-^dy . (14) 
Jo 

In the second integral we transform variables by taking y = 1 — e _l5 ( 1_ *). Here, to is the time at which we have 
exposed the giant component, i.e., when c uncx (t) = c unto (t); since c uncx (t) = 5t 2 and c unex (i) = #te _l5 ( 1_ *), 
to is the smallest positive root of t = e _5 ( 1- *). 
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This integral can be expressed in terms of the exponential integral function Ei(z) [24] and the incomplete 
Gamma function T(a, z), which are defined as 



Ei(* 



-L 



■ dx 



Then, with the integral 
and the Taylor series 



/oo 
x a ^e- x dx . 

J\ acby dy=±(Ei[ac« b ] -Eifae^]) 

^a-*))-^--)-E^T 



taking a = —5(1 — z) and b — — 6 as in (14) gives 

9° hs {z)^ 7 - (Ei[-J(l-«) 



Ei 



-Je-** 1 -* ^!-*) 



°° z m+l 



m=0 

Thus, the coefficients of the observed degree sequence are 



obs _ J_ / 



dx 



to) 



(15) 



(16) 



Now, to approaches e~ s in the limit of large S, and for m < 5, the integral of (16) coincides almost exactly 
with the full Gamma function T(m) since it contains the peak of the integrand. Specifically, in [6] Clauset 
and Moore showed that if m < S — 5 K for some n > 1/2, then 



a 



obs 
m+1 



(l-o(l)) 



r(m) 

5m\ 



1 

Sm 



giving an observed degree sequence of power-law form m 1 up to m ~ 5 and confirming the experimental 
result of Lakhina et al. [14] . 



7 Conclusions 

Having established rigorously that single-source traceroute sampling is biased, thus formally verifying the 
empirical observations of Lahkina et al. [14], and having calculated the precise nature of that bias for a broad 
class of random graphs, there are several natural questions we may now ask. 

Petermann and Dc Los Rios [22] and Clauset and Moore [6] both demonstrated experimentally that when 
the graph does have a power-law degree distribution, traceroute sampling can significantly underestimate 
the exponent a. Although a characterization of this phenomenon is beyond the scope of this paper, it is a 
natural application of our machinery. 

However, a more intriguing question is the following: can we invert Theorem 2, and derive g(z) from 
<7° bs (z)? In other words, can we undo the bias of traceroute sampling, and infer the most likely underlying 
distribution given the observed distribution? Unfortunately, it is not even clear whether the mapping from 
g(z) to g ohs (z) is invertible, and the complexity of our expression for g ohs (z) makes such an inversion appear 
quite difficult. We leave this question for future work. 
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Finally, although several studies claim that using additional sources in mapping the Internet has only a 
small marginal utility [2, 21], Clauset and Moore [6] showed empirically that in power-law random graphs, 
the number of sources required to compensate for the bias in traceroute sampling grows linearly with the 
mean degree of the network. However, a rigorous analysis of multiple sources seems quite difficult, since 
the events that a given edge appears in BFS trees with different roots are highly correlated. We leave the 
generalization of our results to traceroute sampling with multiple sources for future work. 
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